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Abstract. - We derive an expression for the transport mean free path t]_ associated with 
magneto-transverse light diffusion for a random collection of Faraday-active Mie scatterers. 
This expression relates the magneto-transverse diffusion in multiple scattering directly to the 
magneto-transverse scattering of a single scatterer. 



Magneto-transverse light diffusion - more popularly known as the "Photonic Hall Effect" 
(PHE) - has been predicted theoretically some three years ago jij, and has been confirmed 
experimentally one year later [0 . Phenomenologically, the effect has many similarities to the 
well-known electronic Hall effect: given a diffusion current J, the presence of an external and 
constant magnetic field creates a flow in the "magneto-transverse" direction x J, with 
q the charge of the current carriers. By the non-existence of photon charge, the physics of 
the PHE seems different, and compares perhaps better to the so-called Beenakker-Senftleben 
effect in dilute gazes ||. The evident driver behind the electronic Hall effect is the Lorentz 
force acting on a charged particle while colliding with the impurities. According to Ref. [Q 
the PHE finds its origin in the Faraday effect for dielectric scatterers, that slightly changes 
their scattering amplitude. The charge q is replaced by a material parameter V with the 
symmetry of charge, quantifying the Faraday effect of the particle's material. In a homogeneous 
medium the Faraday effect implies a rotation VB per unit length of the polarization vector 
of linearly polarized light. Two other magneto-optical effects in multiple scattering, such as 
the suppression of coherent backscattering in a magnetic field [Q, || , and "Photonic Magneto- 
Resistance" [D are known to exist, and basically originate from the same Faraday effect. 

In an isotropic medium, Fick's phenomcnological law relates the diffusion current to the 
energy-density gradient Vp according to J = — DoVp Qj. -Do is the conventional diffusion 
constant for radiative transfer and is usually related to the transport mean free path £* and 
the transport velocity ve according to D = ^ vgi* . The velocity is relevant only for dynamical 
experiments. Fick's law applies to media much larger than £* and, when supplied by boundary 
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conditions involving the incident flux, can be solved for the emerging current. In a magnetic 
field, the diffusion constant must be replaced by a second-rank tensor. By Onsager's relation 
-Djj(B) = Dji(— B), the part linear in the external magnetic field must be an antisymmetric 
tensor, and Fick's law becomes, 

J = -D(B)-Vp= D Vp D ± B x Vp. (1) 

The term containing D± describes a magneto-transverse diffusion current. In analogy to Do, we 
shall define the transport mean free path £* ± for magneto-transverse diffusion as D± = ^ve£±- 
For the electronic Hall effect, £* ± would be proportional to the Hall conductivity a xy , whose 
sign is determined by the charge of the current carriers. Similarly, £* ± of the PHE can have 
both signs depending on the scattering. 

Anisotropy in the scattering cross-section - quantified by the familiar anisotropy factor 
(cos 0) H - is well known to make the transport mean free path i* different from the extinction 
length £ ||. The latter is the average distance between two subsequent scattering events. In 
this Letter we will show that the PHE can be understood as a generalization of this anisotropy 
factor to magneto light diffusion, which establishes a difference in scattering between "upward" 
and "downward" directions (with respect to the plane of incident light and magnetic field). 
To this end we use our solution for the Faraday-active dielectric sphere |l(| to relate the PHE 
in multiple scattering, quantified by £* ± , directly to the PHE of one single particle. Although 
such a link may be physically clear, it is not evident from previous work Q. A microscopic 
approach provides both the exact sign as well as the role of anisotropy, which will enable us 
to conclude this Letter with a realistic comparison to reported experiments S. 

The magneto-active dielectric sphere has been discussed by Ford et al. JO]. Experiments 
and symmetry arguments show the PHE to be linear in B. Therefore we have derived a linear 
perturbation formula for the scattering amplitude T pp /(B) of one magneto-active dielectric 
sphere [ [To| . This amplitude relates the scattered field in direction p' to the incoming plane wave 
with wave vector p. The differential cross-section, proportional to the modulus squared of the 
scattering amplitude, must satisfy the reciprocity relation da/dfl(p — > p', B) = dcr/d£l(— p' — > 
p, — B). A magneto-cross-section proportional to (p • B) or (p' • B) is parity-forbidden since 
B is a pseudo-vector and p a vector. Together with the rotational symmetry of the sphere it 
must have the form, 

— ^(P - p', B) = F O (0) + det( P , P ', B) Fi(0) , (2) 

0tot dS2 

where cos 9 = p • p' defines the scattering angle 8, cr tot is the total cross-section and 
det(A, B, C) = A • (B x C) the scalar determinant that can be constructed from three vectors. 
The second term in Eq. (g) will be called the magneto cross-section. For a small, Rayleigh 
scatterer one finds F\{&) ~ VBcos9cq/oj. The antisymmetry of this magneto cross-section 
between forward scattering and backscattering causes the PHE to vanish. The Mie solution 
breaks this symmetry and a PHE was seen to emerge J10[ . 

We will use field techniques developed in Refs. |)| |l2fl to calculate the magneto-transverse 
transport mean free path for a random collection of identical Faraday-active dielectric spheres. 
The four-rank tensor Lyjy iPP ' (q) linearly connects field correlations (Ei (u, p + q/2) Ej (uj, p — q 
of incident and outgoing fields in space; q is the Fourier variable of space and p is the optical 
wave vector. On long length scales (q — ► 0) and without absorption it takes the diffuse form, 

T ( n VK\ 1 ik(P,q,B)l lj (~p , ,-q,-B) 

L ijkl , pp , (q, B) - q.D(B)-q ' (3) 

The symmetric form of the numerator is imposed by the reciprocity principle. Rigorous 
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transport theory yields |Tj| , 

l(p, q, B) = i [G(p, B) - G* (p, B)] - iG(p, B) ■ r(p, q, B) • G* (p, B) . (4) 

G(p, B) denotes the Dyson Green's tensor of the ensemble-averaged electric field to be specified 
later; the asterisk denotes hermitean conjugation in polarization space. We left out explicit 
reference to the optical frequency to. 

The tensor r(p,q, B) is linear in q. In real space, the wavenumber q corresponds to the 
gradient in Fick's law (|l|) . The exact relation between T and the diffusion tensor is [[[J |l2) , 

/i2^ 
-^pTrr(p,q,B), (5) 
p 

This shows that only the trace of T comes in. The term Tr<5r linear in B provides the PHE. 
For a low density n of particles with T-matrix T pp / (B), T obeys the Bcthc-Salpctcr equation 

T(p, q,B)=2(p-q) + n^ T pp , (B) • G(p',B) ■ T(p', q, B) • G* (p', B) • T* pp , (B) . (6) 

p' 

We shall solve this equation for the magneto-active Mie particle, up to linear contributions 
in B. r determines the anisotropy in scattering; without magnetic field it can be ascertained 
that r°(p, q) = 2(p • q)/(l — (cos#)), with (cos#) the anisotropy factor in scattering Hj. 

Eq. (||) shows that T is a hermitean tensor, linear in q. We will ignore the longitudinal 
field in Eq. (^J) that turned out to be of higher order. This makes T as well as G effectively 
a hermitean 2x2 matrix. It is convenient to express all such matrices with respect to a 
helicity base, where er(p) = — <j(— p) = ±1 the helicity of a plane wave with wave vector p. 
They must be a linear combination of the identity U and the three Pauli spin matrices <7 Xi y iZ 
[ |i3| . If the z-axis is taken along the wave vector p, the Green's tensor can be written as 
G aa ,(p,B) = G (p)U + Gi(p) (B .p)<7 B , with G (p) = l/[to 2 /c 2 - p 2 + ito/c £} the Dyson 
Green's function in terms of the extinction length £; furthermore G\(p) ~ Gq. 

We separate the magneto-optical term as T(B) = r° + #r(B). Mirror-symmetry imposes 
that T PP /(B) = T_p_ p /(B) so that, by Eq. (|), T(p,q,B) = -r(-p,q,B). Finally, this 
implies 

<Sr (T(T -(p,q,B)=ai(p-qxB)U + <r z a 2 (B • q) + a 3 (B • p)(p • q) 

+ a 4 (B y q x + B x q y ) cr x + (B y q y - B x q x ) a y (7) 

in terms of 4 real-valued coefficients a n to be determined. They correspond to the 4 Stokes 
parameters I, Q, U and V of the diffuse light. The presence of matrices other than the unit 
matrix in Eq. (|t]) implies that the diffuse radiation (|4|) becomes polarized in the presence of a 
magnetic field. By Eq. (||), the PHE is determined by Tr^r ~ aj., for which Eq. (^) provides 
a closed equation, as we will show now. 

The linearization of Eq. (||) in the magnetic field generates three coupled contributions: 
ST = STi + ST 2 + ST 3 . The first comes from ST(B) itself, the second from 5T pp ,(B) and 
the third from <5G(B). The T-matrix of a Mie sphere without magnetic field depends on the 
scattering angle 9 and the azimuthal angle <f> in a frame where the incident wave vector is 
directed along the 2;-axis|§|. With respect to a helicity base, it takes the form, 

^ P „y = -— e-** [Si(0) + S 2 (0)<7(p)<r' (p')] - -—e-^ [(Si + S 2 ) U + (Si - S 2 )a x ] . 

(8) 
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The bars denote complex conjugation and stem from our different convention than the one 
used in Ref. ||. Since a x a y = icr z and u z a x = ia y many terms generated by ST in Eq. 
contain one of the traceless Pauli spin matrices. Since SG ~ <j z it follows that Tr<5r3 = 0. 
Since a\ — U one term involving 04 survives in Tr <5ri but is eliminated by the integral over 
the azimuthal angle <f>. The first contribution simplifies to, 



Tr£ri(p,q,B) = J d 2 p' F (6')Tr <5r(p'. q, B) = 2ai (cos9) (p • q x B) . (9) 

with F (9) defined in Eq. (0). The second equality follows by substituting Tr <5r(p', B, q) = 
2ai(p'-qxB). The contribution Tr<5r 2 is obtained by substituting T° = 2(p-q)/(l — 

(cos 6)), 

Tr^r 2 (p,q,B) = - /d 2 p'(p'.q)det(p, P ',B)F 1 (0), (10) 

1 - (cos V) J 

and involves the magneto cross-section introduced in Eq. (j|). From Tr 5T^ = 0, and Eq. (Q) we 
have Tr^ri = 2ai(p • q x B) — Tr<5r2- Eqs. (^|) and ( |l0| ) hold for any q and it is convenient 
to choose q = p x B. This yields the desired equation, 

ai = Tl / ... 2 [d 2 p> det y>V ,B) tt r d si ^0 Fl{e) , (ii) 

(1- (cosfl)) 2 J | P xB| 2 (1 - (cos6»)) 2 7 

The last equality follows upon integration over the azimuthal angle <j>. 

We can now insert our result for TrT into formula (|5|) for the diffusion tensor. Since 
Tri5r = 2aiPieykqjBk, the final result reads, 

This equation identifies I* = £/(l — (cos 9)) as the transport mean free path of light in 
multiple scattering, and t*> = a\l as the transport mean free path for magneto-transverse 
light scattering. In Ref. jl0| we have shown that the integral in Eq. (|ll| ) equals exactly the 
normalized PHE of one magneto-active Mie sphere. Eq. ( |l2| ) thus states that the PHE in 
multiple scattering is directly proportional to the normalized PHE of one single Mie sphere, 
including the same sign, and amplified by the factor 1/(1 — (cos#)) 2 . Note that the magneto- 
transverse transport mean free depends even more on the anisotropy factor (cos 9) than the 
conventional transport mean free path £*. Unfortunately, we have no simple explanation for 
this. 

In Fig. 1 we show £* ± /£* as a function of the size parameter, for an index of refraction 
rn = 1.128, corresponding to CeF3 in glycerol. Around x « 40, (radius 2 /im) we calculate 
tjt = +0.06 VBX which, for V = -1100 rad/mT (at temperature T = 77 K) and vacuum 
wavelength A = 0.457 /im yields £]_/£* = -2 • 10~ 5 /T. The experimental value is £±/£* ps 
— 1.1 ± 0.3 • 10~ 5 /T for a 10 vol-% suspension We estimate a systematic error of at least a 
factor of two in attributing values to £± and £* on the basis of Fick's law. Another uncertain 
factor is the broad size distribution in this sample, which probably washes out the oscillations 
of £* ± /£* as a function of size parameter x that are seen in Fig. 1. Besides CeF3, the present 
theory is also able to reproduce the measured sign and magnitude for the PHE of polydisperse 
samples containing ZnS, AI2O3, TiC>2 and EuF2. The sign can change as a function of x and 
m and cannot be predicted by a simple argument known to us. The estimated anisotropy 
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Fig. 1. - Ratio of magneto-transverse transport mean free path and isotropic transport mean free 
path £* ± /£* = (1 — (cosf?))ai, as a function of the size parameter x — 2na/\ of the Mie spheres. The 
index of refraction m — 1.128 of the spheres corresponds to CeF3 in glycerol. The ripple structure is 
attributed to resonances in the various partial waves that built up the scattering amplitude for large 
x. For size parameters x < 1 (Rayleigh regime) we find that £±/£* ~ x 5 vanishes rapidly. Dashed 
line: the same for m —* 1 (Rayleigh-Gans regime), for which a finite value is seen to survive. 

factor for our sample equals (cos 6) « 0.9. The amplification factor 1/(1 — (cos6>)) 2 for the 
PHE is therefore significant, and improves theoretical predictions made for Rayleigh scatterers 
considerably g. 

We conclude that the present theory for spheres of arbitrary size is a significant step forward 
in the qualitative and quantitative understanding of the magneto-transverse light diffusion. We 
thank Geert Rikken, Anja Sparenberg and Venkatesh Gopal for the many discussions. This 
work has been made possible by the Groupement de Recherches POAN. 
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